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The measurement of frequency shifts for light beams exchanged between two test masses nearly 
in free fall is at the heart of gravitational wave detection. It is envisaged that the derivative of the 
frequency shift is in fact limited by differential forces acting on those test masses. We calculate 
the derivative of the frequency shift with a fully covariant, gauge-independent and coordinate-free 
method. This method is general and does not require a congruence of nearby beams’ null geodesics 
as done in previous work. We show that the derivative of the parallel transport is the only means by 
which gravitational effects shows up in the frequency shift. This contribution is given as an integral 
of the Riemann tensor -the only physical observable of curvature- along the beam’s geodesic. The 
remaining contributions are: the difference of velocities, the difference of non-gravitational forces, 
and finally fictitious forces, either locally at the test masses or non-locally integrated along the beam’s 
geodesic. As an application relevant to gravitational wave detection, we work out the frequency shift 
in the local Lorentz frame of nearby geodesics. 


Introduction- The exchange of light beams between 
(almost) free falling test masses and the measurement of 
the corresponding frequency shifts (see Fig. 1) is at the 
heart of any thought (real) experiment devised for mea¬ 
suring in principle (in practice) space-time curvature. It 
is in fact the key for the direct observation of gravita¬ 
tional waves (GW) by interferometer detectors on the 
ground [1, 2] and in space [3], and pulsar timing arrays 
[4] . Particularly at low frequency, GW detectors are lim¬ 
ited by differential forces acting on the test masses [5] 
and, as such, the derivative of the frequency shift may be 
a good observable of space-time curvature and, in gen¬ 
eral, a means to separate true gravitational forces from 
spurious effects. Recently, two different approaches [6, 7] 
have put forward a formalism that allows the frequency 
shift to be computed in terms of an integrated measure of 
curvature. This is distinct from earlier attempts that of¬ 
ten relied on simplifying assumptions, e.g. metric expan¬ 
sion, geodesic deviation, choice of a preferred coordinate 
system, or fixing an a priori gauge (see the introduction 
of Ref. [6] and references therein). 

The only physical covariant quantity that unambigu¬ 
ously describes the effect of curvature in vacuum is the 
Riemann tensor [8]. In effect, other general relativistic 
variables, such as the Ricci tensor and the Ricci scalar 
are identically zero in vacuum, even in a curved space- 
time. Additionally, the Christoffel symbols can be set to 
zero by a proper change of reference frame and the met¬ 
ric itself is in general gauge dependent [9] . Consequently, 
all those quantities are not good observables of the true 
gravitational effect and, as such, they might eventually 
lead to ambiguous results. Previous work has already 
pointed out, although in different formulations, that the 
frequency shift is sensitive to an integrated measure of 
the Riemann tensor over the space-time between the two 
test masses. Those formulations worked out this con¬ 
tribution by defining either a null congruence of present 



FIG. 1. Instantaneous Minkowski diagram for the thought experi¬ 
ment of two test masses exchanging light beams and measuring the 
corresponding frequency shift. The 45° dot-dashed line is the null 
geodesic connecting the emission event e to the reception event r. 
Null geodesics at later instants are shown as thinner lines. The 
other dashed lines (with slopes > 45°) are the two time-like emit¬ 
ter and receiver’s geodesics. 4-velocities, 4-forces and the beam’s 
4-momentum are also displayed for clarity. The derivative of the 
fractional frequency shift depends on the difference between ten¬ 
sors, including the Riemann tensor and fictitious forces, at r and 
e, the latter being delayed by the beam’s light time along the null 
geodesic. 

and past null geodesics [6], or a non-standard time-like 
congruence of the emitter’s velocity and the receiver’s 
velocity [7]. Although the relation between those two 
approaches is not fully clear at the moment, the results 
are similar. 

In this work we give a covariant gauge-independent so¬ 
lution to the problem that does not need the definition 
of a congruence of curves. Instead, starting from funda¬ 
mental principles, in particular the parallel transport of 
4-vectors, we derive the gravitational contribution to the 
frequency shift in a very natural way. It is in fact reason¬ 
able to expect that the gravitational effect should show 
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up quite naturally and directly from the physical inter¬ 
pretation of the formulae. In doing so, we split the con¬ 
tribution of gravity from those effects that can be made 
zero under some reasonable assumptions or coordinate 
transformation -the fictitious forces. 

This paper is structured as follows. We calculate the 
derivative of the frequency shift and make the parallel 
transport explicit in this formalism. We examine the 
terms arising from the differentiation in more detail, and 
we show that the gravitational effect comes from the 
derivative of the parallel transport. Finally, we discuss 
the limit for nearby geodesics, interesting for the appli¬ 
cation to GW detection. 

The derivative of the frequency shift- Consider a 
thought experiment (see Fig. 1) where a test mass, the 
emitter, emits light beams toward another test mass, 
the receiver. We shall consider only a single emitted 
beam and do a general relativistic calculation of the mea¬ 
sured frequency shift. The emission (reception) event is 
e = {r = x^). The emitter (receiver) is moving with 
4-velocity (v^) under the action of an external non- 
gravitational 4-force per unit mass [10] f^ (g^) and grav¬ 
ity. The proper time of the emitter (receiver) is Te (r^). 
The light beam is characterised by the 4-momentum 
and the null-geodesic affine parameter A. We want to 
derive information about the underlying space-time cur¬ 
vature from the observed frequency shift. It is worth 


noting that a phasemeter employed in a space-borne GW 
interferometer measures a phase difference between the 
incoming phase and the phase of a local laser beam [11], 
or the light may be split, directed toward different paths 
and then recombined as in a ground interferometer. In 
any case, the measured phase difference is related to the 
frequency difference via Sui = d5(p/dt. Therefore the nat¬ 
ural observable for this calculation is the frequency differ¬ 
ence between the received frequency Wr and the emitted 
frequency Wg: Suj = uj,- — ujg. Now, the emitted (received) 
frequency is oje = fc^(e)u^(e) [wr = fc^(r)v^(r)] [9, 12] 
and their difference can be written as follows 


Soj = fc^(r) [u'^(r) — u^{e) — 6u ^], (1) 


where has been subjected to a parallel transport from 
e to r along the light beam’s null geodesic [13], i.e. 

= 0. Here, compared to Ref. [13], the paral¬ 
lel transport has been made explicit such that u^(r) = 
u^{e) + 5u^, where 


I r^^(A)u“(A)fc'5(A)dA 


( 2 ) 


can be found by integrating the parallel transport equa¬ 
tion. Also, Ar = A(r) and Ae = A(e). Additionally, as 
r and e are causally connected by the beam’s geodesic, 
it turns out that A and Te are in effect functions of r^. 
Differentiating Eq. (1) with respect to Tr, the result is 


d^w 

dTr 


Dfc 


^ K(r) - u^(e) - Ju'"]-bA:^(r) ^'"(r) -/'"(e) - fc^(r 
aTr L J 


DSu^ 


diff. velocity 
+ par. transport 
+ fictitious 
(a) 


■v 

diff. force 
+ fictitious 

(b) 


par. transport deriv. 
+ fictitious 

(c) 


( 3 ) 


Working out the differential-velocity term (a)- This 
term accounts for the difference of velocities, u'"(r) — 
M'"(e), but also the parallel transport both projected 
along Dfc^/dTi-. In the local Lorentz frame (LLF) at r, 
= 0 so Dfc^/dTr = dk^/dr,. but, being an integral 
along a finite curve, 6u^ d. If the emitter and re¬ 
ceiver are assumed to be infinitely close to each other 
such that the integral can be evaluated at a common ref¬ 
erence frame coincident to the above-defined LLF at r 
(we shall later call it the LLF for nearby geodesics) then 
additionally 6u^ = 0. Also, for every point along the 
null geodesic an LLF can also be defined such that the 
integrand in Eq. (2) will be zero locally. Therefore, while 
Dfc^/drr includes local fictitious forces at r, instead Su^ 
reproduces the non-local effect coming from the integra¬ 
tion of an infinitesimal fictitious force along the beam’s 
null geodesic. 


Working out the differential-force term (b)- Two con¬ 
tributions are included, the real force on the receiver, 
g^{r), and an effective force on the emitter, /'"(e), which 
is given by 

/' = I/'+ r;x(»^(4) 

where dre/dr^ = 1 -\- Sui/uje. Because the non- 

gravitational forces are assumed small, it is safe to ap¬ 
proximate dre/dr^ ~ 1 to first order. The effective force 
/'"(e) contains a term evaluated at e that, in general, 
is not zero when only the receiver’s LLF is considered. 
However, this becomes zero in the nearby-geodesics LLF. 
Therefore this second term may be interpreted as a local 
fictitious force at the emitter. 

Working out the parallel transport derivative term (c)- 
This is the gravitational contribution that the deriva- 
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tive of the frequency shift is sensitive to. Before going 
through the calculation, it is worth noting that may 
be still treated as a second rank tensor when a is fixed 
[14], i.e. Fg^ transforms as a tensor for fixed a = d, 

under a general transformation 
from the unprimed to the primed coordinate system. In 
differentiating the integral in 5u^ with respect to Tj-, the 
extremes of integration depend on r^, so the commutation 
between the derivative and the integral gives 

e 

( 5 ) 

where the integrand depends implicitly on A. The last 
contribution becomes zero in the nearby-geodesics LLF, 
thus it acts as local inertial forces at e and r. From now 
on, the focus is on rearranging the integral in a more 
familiar form. Applying the covariant derivative, the in¬ 
tegrand may be recast as follows 

/ N (6) 

= + rf.ri J + F^^d. {u^^kP) vr 


drv 


r-Ar 


dA 


[KpU-kP) dA+ 


the same. So, looking at Eq. (3), the role of the emitter 
and the receiver may be swapped, and in particular in 
the integral, v'^ O and roe. The integral will for¬ 
mally be the same, but with reversed integration limits. 
Because of this symmetry, the original derivative of the 
parallel transport can be rewritten as 


_ 1 

dr,. 2 ,/x 




+ dA. 


( 7 ) 


where the Riemann tensor arises from the antisymme¬ 
try between a and v. The second term can be recast 
as \yykP — u'^v°‘) + kPCvU'^] , where = 

v''k/„u°‘ — u''Vi,v°‘ = —£u'y“- It turns out that: (i) it 
is zero in the LLF along the null geodesic, (m) it is zero 
when the two velocities are equal, and (in) it is zero 
in the nearby-geodesics LLF. Integrated out along the 
beam’s geodesic, it therefore gives a non-local fictitious 
force. Instead, ^ d can not be set to zero in a 

purposely chosen LLF: it is in fact an intrinsic property 
of space-time that quantifies curvature. This curvature 
contribution is 


In effect, the quantity between parentheses reproduces 
the first two (positive) terms of the Riemann tensor. To 
make it more explicit, it is useful to consider a thought 
experiment where the role of the emitter and the receiver 
are reversed. The observed physical effect will be exactly 


n^ = ll\'^,^pu‘^vPk’^dX. ( 8 ) 

Under the general assumption that the non-gravitational 
forces are small, then the original formula becomes 


d^a; 

dr. 


dTr 


[u'^(r) 


-u^(e)] +kf,{r) [g^^{r) 


rie)] + k^ir)n^ + 7flct [r(r), F(e), F(A)], 


( 9 ) 


where qgct collects all additional terms that may depend 
on F[(^ either locally (at event e or r) or non-locally (via 
an integral along the beam’s null geodesic). In any case, 
7fict goes to zero in the nearby-geodesics LLF. The above 
result shows that not only the derivative of the frequency 
shift depends on the relative velocity and relative non- 
gravitational acceleration, but also contains a contribu¬ 
tion from the integral of the Riemann tensor along the 
beam’s geodesic, plus local and non-local fictitious forces. 

The local reference frame for nearby geodesics- As an 
application, it is worth calculating the limit for nearby 
geodesics where a local approximately-inertial reference 
frame common to both the receiver and the emitter can 
be set up. In such conditions, dr, ~ dre ~ dt so the 
effect of evaluating, for instance, a 4-vector at emission 
or reception is just a light time delay St = Sx/c, where 
Sx is the nominal distance between the emitter and the 
receiver. Additionally, k^ = lo^jcn^^ 5x^ = Sxn^ and 
= (I,h). The Minkowski diagram of Fig. 1 shows 
this LLF and the geodesics of the emitter, receiver and 


light beam. In this reference frame the receiver appears 
still, and if we further assume that the emitter’s apparent 
velocity is negligible to first order, then the derivative of 
the fractional frequency shift becomes 


d 

dt \iOeJ 


= • [v{t) - uit - St)] 

+ ■ Wit) - fit - <^ 0 ] 

+ ^Sx [Rfj,uoeit) - Rf,vmit - <5t)] u'" n'', 


( 10 ) 


where dh/dt = ft x h and U is the observed rotation 
of the line of sight [7]. The additional fictitious forces 
contained in the more general result of Eq. (9) are in¬ 
tentionally neglected here. As per the 4-velocity and 4- 
forces, the Riemann tensor must also be evaluated as a 
difference at the reception time and the delayed emission 
time. Moreover, all tensors are here evaluated in the de¬ 
fined LLF, and so is, of course, the Riemann tensor. So, if 
we want to relate this Riemann tensor to the one defined 
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in the so-called transverse-traceless gauge wave frame of 
an incoming GW, RJqjq = 1/2/iJj^Qg, a rotation will 
naturally show up transforming the wave frame into the 
LLF. 

Conclusions- This paper has addressed the issue of in¬ 
timately understanding how space-time curvature can af¬ 
fect the frequency shift of light beams exchanged between 
two test masses. The nature of the problem is general, 
but ultimately relevant to GW astronomy. Gompared 
to previous work, the calculation is straightforward and 
does not require a congruence of curves, but is always 
matched to the ultimate physical interpretation. 

Without assuming any coordinate system or fixing the 
gauge, thus retaining the full covariant nature of the for¬ 
malism, we have found that the only gravitational contri¬ 
bution to the time derivative of the observed frequency 
shift comes just from the derivative of the parallel trans¬ 
port of the emitter’s 4-velocity along the light beam’s 
geodesic to the receiver. This is, in turn, given by an 
integral of the Riemann tensor along the null geodesic, 
although in a different form compared to previous work 
[6, 7]. Additionally, fictitious forces show up locally at the 
receiver, at the emitter and non-locally integrated along 
the beam’s geodesic. Those terms were not considered 
before, although we note here that they are all propor¬ 
tional to thus they are negligible under reasonable 
assumptions that are always met in practice. 

As an application, we have calculated the derivative 
of the frequency shift for nearby geodesics in the refer¬ 
ence frame of the receiver, showing that it is given in 
terms of differences between velocities, non-gravitational 
forces, and the Riemann tensor, all evaluated at time of 
reception and delayed time of emission. Additionally, the 
observed Riemann tensor can be directly related to the 
usual Riemann tensor defined in the wave frame where 
the gauge is transverse traceless. 

It is worth commenting that all formulae containing 
the metric perturbation to an underlying flat metric, 
as customary in GW astronomy (see e.g. Ref. [15]), can 
be derived from the result contained in this paper by 
choosing a convenient coordinate system and by fixing 
the gauge. Nonetheless this does not alleviate, at least in 
principle, the problem of fictitious forces that may still be 
dependent on itself. But, ultimately, the existence of 
a local LLF for both the receiver and the emitter justifies 


the use of time-delayed differences in the calculation of 
the detector response to GWs. An example of which is 
found in Ref. [5] where the detector response was given in 
terms of acceleration, a quantity that is very convenient 
for marginalising over detector systematics. 
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